Effect of Electron-Electron Interactions on Rashba-like and Spin- Split Systems 
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The role of electron-electron interactions is analyzed for Rashba-like and spin-split systems within 
a tight-binding single-band Hubbard model with on-site and all nearest-neighbor matrix elements 
of the Coulomb interaction. By Rashba-like systems we refer to the Dresselhaus and Rashba spin- 
orbit coupled phases, while spin-split systems have spin-up and spin-down Fermi surfaces shifted 
relative to each other. Both systems break parity but preserve time-reversal symmetry. They 
belong to a class of symmetry-breaking ground states that satisfy: (i) electron crystal momentum is 
a good quantum number (ii) these states have no net magnetic moment and (iii) their distribution 
of 'polarized spin' in momentum space breaks the lattice symmetry. For all members of this class, 
the relevant Coulomb matrix elements are found to be nearest-neighbor exchange J, pair-hopping J' 
and nearest-neighbor repulsion V. These ground states lower their energy most effectively through 
J, hence we name them Class J states. The competing effects of V — J' on the direct and exchange 
energies determine the relative stability of Class J states. We show that the spin-split and Rashba- 
like phases are the most favored ground states within Class J because they have the minimum 
anisotropy in 'polarized spin'. We analyze these two states on a square lattice and find that the 
spin-split phase is always favored for near-empty bands; above a critical filling, we predict a transition 
from the paramagnetic to the Rashba-like phase at a critical J(J c i) and a second transition from 
the Rashba-like to the spin-split state at J C 2 > J c i- An energetic comparison with ferromagnetism 
highlights the importance of the role of V in the stability of Class J states. We discuss the relevance 
of our results to (i) the a and /3 phases proposed by Wu and Zhang in the Fermi Liquid formalism 
and (ii) experimental observations of spin-orbit splitting in Au(lll) surface states. 

PACS numbers: 



I. INTRODUCTION 

Numerous broken-symmetry phases have been pre- 
dicted to arise from electron-electron interactions. These 
phases are often framed in the language of Pomeranchuk 
instabilities of the Fermi liquid pJiiJ. In this language, 
the interactions that drive a symmetry breaking are 
purely phenomenological, though work has been done 
to derive the Landau interaction parameters from a 
more physical model, e.g. the Hubbard model; in these 
works [5|, l6| particular attention has been drawn to the 
Hubbard parameter U (describing on-site repulsion) 
because of the strong interest in high-temperature 
d-wave superconductivity There has been further 
progress made by Lamas and coworkers in understanding 
how the lattice affects FL instabilities in two dimensions 
Q, focusing on Hubbard interactions U and V (nearest- 
neighbor repulsion). 

In the tight-binding representation, U and V are 
matrix elements of the Coulomb interaction that are 
diagonal in the density operators; there are also non- 
diagonal matrix elements |9Hl3j. In total, there are four 
nearest-neighbor interactions, denoted by V (diagonal) 
and J, J' and At (off-diagonal). This Hamiltonian 
(including U) is usually termed 'generalized Hubbard 
model'. 

In this paper, we ask: other than ferromagnetism and 
antiferromagnetism, what kinds of broken-symmetry 
phases involving spin arise in the generalized Hubbard 



model? It was initially proposed by one of us [14 1 
that the interaction J (nearest-neighbor ferromagnetic 
exchange) plays a crucial role in displacing the Fermi 
surfaces of spin up and down electrons relative to one 
another; the resultant spin-split metallic phase is illus- 
trated in Fig. (1-a). Subsequently, the spin-split phase 
has also been predicted to arise in the FL framework by 
Wu and coworkers and by Chubukov and Maslov 
filij . A possible realization of the spin-split phase was 
originally proposed to be the low-temperature broken- 
symmetry phase of Chromium [3]; recently Varma and 
Zhu have proposed that it may describe the 'hidden or- 
der' phase in the heavy fermion compound U RU2S12 [l7| . 

In addition to the spin-split phase (dubbed the a 
phase by Wu et al), Wu et al also showed that many 
spin-orbit coupled phases (the /3 phase) may arise from 
a phenomelogical Landau interaction in the p-wave spin 
channel. In two dimensions, these phases include the 
Rashba [l8j], Dresselhaus [l9| and helicity spin-orbit 
coupled systems; for convenience, we call these group 
of spin-orbit coupled phases Rashba-like. The Rashba- 
polarized state is illustrated in Fig. (1-b). 

In this paper we return to the tight-binding Hubbard 
formalism and analyze the effects of on-site repulsion 
and all nearest-neighbor matrix elements of the Coulomb 
interaction. We demonstrate that the spin-split state 
and the Rashba-like phases are part of a broader class of 
symmetry-breaking ground states for which the relevant 
matrix elements are nearest-neighbor exchange J, pair 
hopping J' and nearest-neighbor repulsion V. This 
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FIG. 1: (a) Fermi energy contours of a spin-split state (b) 
Fermi energy contours of a Rashba-polarized state 



class of ground states has three properties: (i) electron 
crystal momentum is a good quantum number; (ii) 
these states have no net magnetic moment; (iii) their 
distribution of 'polarized spin' in momentum space 
breaks the lattice symmetry. In this class of ground 
states, we compare the stability of different phases 
and examine the individual roles played by J, J 1 and 
V in changing the dispersion and the ground state energy. 

In our analysis of the Rashba-polarized state, we are 
also motivated by the following observation. Photoc- 
mission experiments on Au(lll) surfaces have revealed 
a splitting of Rashba-polarized Fermi surfaces; the ex- 
perimental results may be reproduced in a free-electron 



model by a phenomelogical interaction 

(n x k) 



H, 



a R cr ■ 



(1) 



with an electric field perpendicular to the surface (par- 
allel to n) [20U241 I . The Fermi wavevector in Au is k F ~ 
0.16 A -1 and the measured energy splitting of the Fermi 
surfaces is ~ 0.1 leV, corresponding to a Rashba cou- 
pling constant a R ~ Q.34eVA _I . The Rashba coupling 
constant derived from the Dirac Hamiltonian in the non- 
rclativistic limit is 



eft 2 



E 



(2) 



with e and m e the electron charge and mass and E the 
electric field normal to the surface. For the experimen- 
tally observed splitting, the magnitude of the electric field 
in Eq. J2) is 



-eE = 91,228eV/A" 



(3) 



Such a large value of the electric field is justified only 
with the assumption that the conduction electrons 
get very close to the Au nucleus; this is incompatible 
with a free-electron model. In this paper we ask the 
question: can electron-electron interactions cause a large 
enhancement of the single-particle spin-orbit splitting so 
that the observed Fermi surface splitting is explained 
with a much weaker electric field than the value in Eq. 
©? 

This paper is organized as follows. In Sec. II, we 
discuss the relevance of various Hubbard parameters and 
identify a class of ground states for which the relevant 
parameters are J, J' and V. With these important 
parameters, we derive a reduced Hamiltonian. We solve 
this reduced Hamiltonian for Rashba-like phases in Sec. 
Ill and for the spin-split state in Sec. IV. In Sec. V, we 
find that V plays a pivotal role in the phase transition 
between the Rashba-like and the spin-split state. In 
Sec. VI, we consider the stability of Class J ground 
states other than the spin-split and Rashba-like phases; 
we also discuss various symmetries of the class. In Sec. 
VII, we investigate the stability of the spin-split and 
Rashba-like states in comparison with ferromagnetism. 
Sec. VIII draws the connection between the formalism 
in this paper and the Fermi liquid a and /3 phases [15|. 
In Sec. IX, we summarize the key physical ideas and 
discuss possible generalizations to other lattices. Several 
technical points are discussed in the Appendix. 



II. REDUCED HAMILTONIAN 

We consider a single band generalized Hubbard model 
with on-site and all nearest- neighbor interactions [25| : 



H = H t +H v + Hj + Hj 



(4) 
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Respectively, the individual terms are the nearest- 
neighbor hopping and the various electron-electron 
interactions that correspond to on-site repulsion U, 
nearest- neighbor exchange J, pair- hopping J', nearest- 
neighbor repulsion V and correlated hopping At. All the 
electron-electron interaction parameters are positive [l2|. 

We consider a two-dimensional square lattice and as- 
sume that the spinor is confined in the x-y plane; we 
denote this restricted set of bases by the index {f2 = ±1} 



C ka =-j= [ ' k] + { - } - 'I/.:! <k , 



(5) 



where we choose the convention that the basis {a =J [,i} 
is quantized in the z direction. We write the momentum- 
dependent amplitude ip as ((k) + ix{k), with ((k) and 
x(fe) real. The physical interpretations of £ and x are 
clear - a spinor in the x-y plane may be expressed as 



^2 \{<rx) + * K, 



(6) 



The normalization \ip\ 2 = 1 is just the condition that 
1. 



(a x ) 2 + (a v ) 2 



In assuming that the spin orientations of all electrons 
are coplanar, we are anticipating an exchange interaction 
that favors parallel alignment of spins. An out-of-plane 
degree of freedom in spin space will tend to disrupt 
the parallel alignment of spins and is energetically 
suppressed. 

In the coplanar basis characterized by ip, the kinetic 
hopping term has the form 



H t = -2 1 y^(cos k x + cos k y ) c ka c kn 



k.n 



(7) 



just like in the a basis. Before transforming to the f2 ba- 
sis, the Coulomb interaction has terms quartic in creation 
and annihilation operators in the a basis: 

C k+q,a C k> -q,-a C k' -oCk,o 

k.k' ,q,S 
a 

~ ON 2—1 C k+q,<j C k'-q,<7< C k' ,<r'C-fc,<T 



J' 



k,k' ,q,8 
a. a' 



~2N 2—1 C k+q,a C k'-q-c 



V 



C-k' , — crCk,<T 



k.k ,q,S 



HV= 2N eiqSc k+q,<r C k'- q ,a> C k',v> C k,<> 
k,k' ,q,8 

H At = ^L J2 {e lkS + e lk ' s + e^ k -^ s + e^'-^H) 

k,k' ,q,8 
a 

X C k+q,o- C k' -q,-<J C k' -a Ch - a ' 



(8) 



We convert this to the il basis and decouple in the 
(cL n Cfe n ) direct channel. Define mean fields 

n k =n k + + n k - = {c\ + c k+ ) + (c f k _c k -) and (9) 

A k =n k+ - = (c k+ c k+ ) - (c^_Cfc_). (10) 

In this paper we consider the set of ground states sat- 
isfying the symmetry constraints 



^2ip(k)A k = and 

k 



j~y fc< V<(fc)A fe = o. 



(ii) 

(12) 



k.8 



We will call them 'Class J ground states'. We deduce 
from Eqs. ([5]) and (fTUj) that Eq. (ITTt sums over the 
spins of all the polarized electrons, hence these ground 
states have no net magnetic moment - ferromagnetism 
does not belong in this class. Eq. (fT2|) ensures that the 
momentum distribution of the 'polarized spin' ip A does 
not have the symmetry of the underlying lattice; in two 
dimensions this symmetry is the point group C nu for a 
lattice with n-fold symmetry. From the choice of mean 
fields in Eqs. © and (ITU)) , we deduce that electron 
crystal momentum k is a good quantum number in Class 
J, hence antiferromagnetic and spin-density-wave states 
are excluded. 

We prove in Sec. VI-C that these ground states 
break parity. For the special case that £(fe) x(&) are 
invariant under k — > — fe, we also show that they preserve 
time-reversal symmetry. 

There are many variations of ground states that fall 
into this class. We consider the following cases: 

(i) A breaks the lattice symmetry but ip does not. 
The simplest example is when ip is independent of 
momentum - the spin orientation of all electrons point 
along the same axis - and there is a relative displacement 
of spin up and down Fermi surfaces with respect to 
each other. This is the spin-split state shown in Fig (1-a). 

(ii) ip breaks the lattice symmetry but A does not. 
The simplest example is an expansion of one Fermi 
surface with respect to the other in a basis that rotates 
spin by 2tt rad as we go around the Fermi contour, 
i.e. winding number of 1. These include eigenstates of 
helicity (spin parallel or antiparallel to momentum) and 
eigenstates of the Rashba and Dresselhaus spin-orbit 
interactions. The Rashba state is shown in Fig (1-b). 

(hi) It is also possible that ip and A both break the 
lattice symmetry. These exotic ground states are shown 
to be energetically unfavorable in Sec. VI-B. 

These constraint will have important implications on 
which Hubbard parameters are relevant to the ener- 
getics of polarization. In particular, we will find that 
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the nearest-neighbor exchange J is the greatest driv- 
ing force for symmetry breaking, hence the name Class J. 



We proceed to reduce the Hamiltonian in Eq. @ . 
introduce matrices 



We 




(13) 



(CC' + XX') 

HxC-Cx') 



i(xC-Cx') 
-(CC' + xx' 



(14) 

To simplify notation, we have written as £ and C(&') 
as C' in Eq- (p~4|) . and similarly for \- For the rest of 
this paper, we use the convention that II = IT(fc') and 
T = T(fc,fc') as defined in Eqs. ^ and CD]}. The 
diagonal and off-diagonal elements of T are proportional 



to ei- 



and 



,er x 



(er') respectively. 



There are two values of momentum transfer q that 
contribute to the mean field: q = and q = k' — k. 
For example, consider the on-site repulsion U . We may 
organize 

J_£{q = k' -k) 



Hu - H$ ° } -r <<, 



After the mean-field decoupling 
U 



tt{i=°) 



(q = fc'-fc) 



(4 

27V £ 



~2N ^ 

k,k\S 

u 



J 

-fc- 



-k- 



k,k',S 



IT 



T 



Cfc+ 
Cfc- 

( c k+ 
\Ck- 



(15) 

(16) 
(17) 



We have neglected constants proportional to (c^c)(c ! c) 
that will affect the expectation value of the ground state 
energy but not the single-particle energy spectrum. Wc 
carry out the integral over k! . Defining n as the total 
number of electrons, ^ fe , II = n 5 rin i which is a constant 
matrix independent of polarization. Eq. (|lll) leads to 
J2k' — 0- We conclude that Hjj does not contribute 
to Class J-type polarization. For the same reasons, we 
r(«=*'-*0 and ff<*=°>. 



may discard H { J 



We also consider 

rr(q= 
11 At 

At 



TT(q = k -fe) 

11 At 



2N 



E ■ 

h,h',8 



Ah- 8 



(4 + 4.) («) ■ 



(18) 

The sum over fc' vanishes because of both symmetry 
constraints in Eq. (fTTj) and (fT2")l . 

The following terms are left: 



Jj(Q=0) 



-y 

2N ^ 



i(k-k')-S ( j 



k + 



k.k'-S 



2n 



Cfc+ 
Cfc- 



(19) 



2iV 



= •' v P *(*+*')-« 



k.k'.S 



(4 



fc+ l 'fe- 



T + n 



Cfc+ 
Cfe- 



(20) 



TT^q—k' — k) 

n v 
V 
2N 



= J_ V f_t 



fc,fc',<5 



T-n] (2D 



11 At 



At 

'2N 



k.k',8 



(e lkS 



ih -8\ 



^k- 



2n 



( c fe+ 

\Ck- 



(22) 



With the exception of iJ^ 4 =0) , the other interactions are 
particle-hole symmetric - we will show that their effects 
are extremized at half-filling. 

We briefly describe the roles of the various Hubbard 
parameters: 

(i) The expectation value of Eq. (|19l) is 



(H ( r 0) ) 



- T 

O.N i—i 



e^-fc)" 5 n k n w . (23) 



k,k' ,S 



The amplitude Yls e lt - k ~ k is negative when ki ~ 
and k[ ±7r or vice versa, hence J favors a separation 
of electrons in momentum space, i.e. polarization . This 
was first pointed out by one of us [26] in connection with 
ferromagnetism. Since Eq. (|23p is independent of the 
choice of basis, we find that J promotes polarization in 
any arbitrary basis; this is true even for ground states 
that do not satisfy the symmetry constraints in Eq. 
(fTT]) and (fi"2j) . Class J ground states are special be- 
cause they reduce their energy most effectively through 
nearest-neighbor exchange J when they break symmetry. 

(ii) The expectation value of Eq. (|2T|) is 



(#<r fc '- fe) ) 



v 

AN 



k,k',S 



,i(h-k'ysl A A 



(24) 



The second term has the opposite effect of J and tends 
to suppress polarization. The first term causes an 
energy splitting between the — ±1 bands and tends to 
lower the energy if polarized electrons with fc « fc have 
parallel spin. This is analogous to the exchange energy 
J that favors parallel spins for localized electrons; 
this was proposed by Heisenberg as a mechanism to 
explain ferromagnetism 27]. The difference is that in 
the momentum representation, J influences the direct 
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energy while V mixes both direct and exchange energy. 

While J unambiguously assists polarization in any 
basis, V affects both direct and exchange energies in a 
manner that is self-cancelling; its net effect on polariza- 
tion depends sensitively on (i) the choice of ground state 
and (ii) band filling. Hence, V plays the crucial role of 
deciding which ground state is energetically favored for 
a given band filling. 

(iii) We show in Sec. III-E that the main effects of J' 
and At are to cancel V and t respectively. 

We now analyze the two simplest examples of Class J 
with this reduced Hamiltonian - they are the Rashba- 
type and the spin-split states. In Sec. VI-B, we 
demonstrate that they are also the most stable in this 
class. 



III. CASE STUDY OF RASHBA-TYPE 
GROUND STATE 

We examine a set of ground states with spin struc- 
tures that diagonalize Hamiltonians of the form H oc 
a x sin k a ± <r y sin fc;,. For specificity, we first consider the 
Rashba spin-orbit basis. For completeness of analysis, 
we add a Rashba single-particle interaction Hr to the 
Hamiltonian in Eq. this originates from the Dirac 
spin-orbit interaction (Eq. ([1])) with an electric field per- 
pendicular to the plane 

H = H t + H R + Hu + Hj + Hj, + H v + H At . (25) 

With only nearest-neighbor hopping, the tight-binding 
version of the single-particle Rashba Hamiltonian (Eq. 
©) is (H 



We introduce the generalized dot and cross products 
sin h 



ki — 



k • k f — ^ kikjSij, 

ij 

fc x In — ^ ^ k j k j €. j j . 



(30) 

(31) 
(32) 



Note that is the Levi-Cevita symbol. The matrices 
introduced in Eq. (fTB"]) and (fT4"]l take the form 



n 



T=A fc 




~(k-k') 
(k x k') 




(33) 
(34) 



We may easily generalize the above results to include 
the tight-binding generalizations of the Dresselhaus and 
helicity spin-orbit bases. The Dresselhaus basis follows 
from changing the spin orientation fk — > sin k x + i sin k y 
in Eqs. (12511 and J2HJ); the helicity basis changes 
fk — > smk x — i sin k y . Without any single-particle spin- 
orbit coupling, all three ground states are degenerate 
because electron-electron interactions cannot distinguish 
between their spin structures. Other than the spin 
structure encoded in all the results that follow in the 
rest of Sec. Ill are applicable to the Dresselhaus and 
helicity states if cxr = 0. 



We first analyze a further reduced Hamiltonian with 
only H t ,H R ,H ( f =0) and H$ =k '- k) in one and two 
dimensions; we defer a discussion of the effects of H l ^ 0} 
and Hj, to Sec. III-E. 



H R =-a R 4 CT ( sin k x (a y ) - sin k y (a x ) I c ka > 

fcerer' V 7 

(26) 

where spin indices a are chosen to be diagonal in the 
z direction; Eq. (J2U) is formally derived in App. A. 
Diagonalizing this Hamiltonian in momentum space, we 
obtain 



H R = a R T)J sm 2 k x +sm 2 k y c^Ckr,- (27) 

k,q 

We distinguish between indices a and rj; the latter de- 
notes the Rashba polarization and takes on values ±1. 
The eigenvectors are 



fk = sin k y + i sin k x . 



(28) 
(29) 



A. One Dimension 

The reduced Hamiltonian is exactly solved in one di- 
mension. Many of the qualitative conclusions generalize 
to two dimensions. 

The generalized dot product reduces to k ■ k' = 
sgn(fe) sgn(fc') and the cross product vanishes. The in- 
teraction part of the Hamiltonian simplifies to 



Hv > j = jj cos ( fc ~ k 'y\r, c ^ 



x (2 J - V)n k , - riVA k , sgn(k) sgn(fc') (35) 



We parametrize the polarization by a dimensionless 
quantity e: n k ± — 1 for fc £ {~k F ± e, +k F ± e}. The 
single-particle energy is 
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E v (k) = — 2 (t — —(2 J — V) sinfc F cose^ cos fc 

+ r) ( a R + — Vsink F sine J sin (36) 



The ground state expectation value of the full Hamil- 
tonian is 



iV 



sin k F a R sin e + 2 1 cos e 



7T 

+ A sin" /,-,, ( 2J cos- < • 1" 



(37) 



From Eqs. (f55|) and (|3"T1) . we see explicitly that J 
promotes polarization by increasing the bandwidth of 
the single-particle energy spectrum. On polarization, 
V simultaneously reduces the bandwidth and increases 
the energy splitting of the Rashba bands; the net 
effect is that the V term in the ground state energy is 
independent of e. 

Minimizing (H) with respect to e, the ground state 
polarization e* is determined by 



a R cose* H — Jsink F sin (2e*) = 2£sine*. (38) 

7T 

This same equation can be obtained by enforcing the self- 
consistency condition 



E„(k F +e*) =E + (k F -e*). 



(39) 



In Sec. III-B, we prove that they are equivalent state- 
ments by minimization with Lagrange multipliers. 

If a R = 0, the polarization is determined by 



t 



cose = — - 



2 J sin k F 

and there is no nonzero solution for J < J c with 
J * 1 

■> C 



2 sin. k F 



(40) 



(41) 



At J > J c there is a symmetry breaking: e may take 
on values ±e*. Note that the effect of J is maximized 
at half-filling; this is consistent with the particle-hole 
symmetry of the nearest-neighbor exchange interaction. 

Up to an arbitrary spin rotation, the Rashba state is 
equivalent to the spin-split state in one dimension [lij : 
we show this in Sec. IV- A. 
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(b) 



FIG. 2: (a) Fermi energy contours of the two Rashba bands. 
Arrows refer to spin directions, (b) Dispersion at fixed k y = 
0. Parameters: t = 1, V = 2, J = 3,a R — 0. This band is 
half-filled. Symmetry breaking occurs for J > J c = 2.5. 



B. Two Dimensions 

The off-diagonal terms in Hy.j are proportional to the 
integral 



j-y(h-«0-« ^ x Afc , 



(42) 



S,k' 



which vanishes if the splitting satisfies = A_^ and 
A/- k — A±k ±k x - The latter condition reflects the 
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symmetry of the underlying square lattice. 
We define the geometric factors 



a =jr y^fcos k x + cos k y ) n k , (43) 

k 

/? =4 X! V s in2 k x + sin 2 k y A k . (44) 



Up to a negative constant — 2t, a is the sum of the 
kinetic energies of all the electrons; j3 can be interpreted 
as the order parameter and is defined to be negative for 
positive a R . We also define a as the value of a when 
there is zero polarization; /3 D = 0. 

The single-particle energy spectrum is then 
1 



E v (k) 



2 | / - -(2 J - V)a ) (cosfca; + cosk y ) 



+ 77 ( a R - -VP ) y 1 sin 2 k x + sin 2 k y (45) 



We observe that V results in an effective momentum- 
dependent Zeeman field that is maximized at momenta 
(±7r/2, ±7r/2) while f3 is the overlap of the polarization 
Afc with this Zeeman field. Consequently, V reduces the 
exchange energy most effectively at half-filling, where 
any polarization takes maximal advantage of the Zeeman 
field. 



The ground state energy is 



N 



-2ta 



a R P + ija 2 



/(a 2 + /3 2 ). (46) 



a {(3) increases (decreases) as the system polarizes; by 
inspection of Eq. (I46[) , we confirm that J and a R favor 
polarization. Whether V assists polarization or not 
is not clear from Eq. (|46l) alone. In one dimension, 
a 2 + p 2 is a constant that is independent of polarization 



(c.f. Eq. (|37p). hence V will not tip the balance in 
either direction. In Sec. III-C and V, we investigate the 
behavior of a 2 + p 2 in two dimensions. 

We wish to minimize the ground state energy with the 
constraint of fixed particle number n. Hence we define 



g=^-m(i>-») 



and impose the condition 
dG 



dn 



= 0. 



(47) 



(48) 



Since we are making an infinitesimal variation to the 
Fermi surfaces of both Rashba bands, Eq. (|48[) only ap- 
plies to states at the two Fermi surfaces. From Eq. 



(E v (k)) 



FS,rj 



(49) 



or 

(E + (k)) FS+ = (E_(k')) FS _ 
which is a generalization of Eq. (|39l) . 
We define a and b: 



(50) 



Ekr) = — a(cos k x + cos k y ) + r\ by sin k x + sin k y . (51) 
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(b) 



FIG. 3: (a) Fermi energy contours of the two Rashba bands 
for J = (b) Fermi energy contours of the two Rashba bands 
for J = 2. Filling = 0.5. Parameters: i = 1, V = 0, a B = 0.2. 
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b is the energy splitting and a is a quarter of the 
bandwidth; for convenience, we call a the bandwidth. 

Given parameters t, J, V, a R and a filling n, Eqs. 
([§]). ([l % (|i3 ]l . (|13]) . (@1J and j^} may be solved numer- 
ically. A few representative results are presented: (i) in 
Fig. (2-a) and (2-b), a symmetry-breaking solution exists 
for a R = and J > J c (ii) in Fig. (3-a) and (3-b), a com- 
parison is made between two Rashba-polarized systems 
with equal a R but different J. 



Eq. (1551) is a good approximation even when extended 
to systems that are medium-polarized (< 50% of 
the total filling); this is illustrated in Fig. (4). In 
Sec. V-A, we provide a geometric explanation of why 
this approximation is especially good near quarter-filling. 

Employing Eq. (|55|) . we cast Eq. (|54"|) in the useful 
form 



I 2 

J = sOi R + —t. 
p a 



(56) 



C. A useful analytic approximation 

We are motivated by the simplicity of the equations 
in one dimension to look for a suitable approximation 
of the two-dimensional equations; this approximation 
permits us to extract an analytic formula that relates 
J, a R and t. We numerically explore the limits of validity 
of this approximation. Where the approximation is not 
valid, we believe our conclusions are at least qualitatively 
correct. 

In one dimension, the Rashba bands can polarize in a 
single direction. Hence, we minimize (H) with respect 
to a single polarization parameter e to obtain the unique 
ground state with polarization e*. In addition, 



-sin k F cos e 



-sin k F sine 



(52) 
(53) 



and the combination a 2 1D + j3f D is independent of 
polarization parameter e* and only dependent on the 
the number of electrons; this is a consequence of number 
conservation. 

In two dimensions, the Rashba bands can polarize 
anisotropically, hence we have to minimize (H) with re- 
spect to all the occupation numbers on the Rashba Fermi 
surfaces (see Eq. (|48p). This minimization process is 
equivalent to minimizing (H) with respect to mean fields 
a and /3; this claim is justified by the Hohenberg-Kohn 
theorem [29| in App. B. After minimizing Eq. ((46]), we 
obtain 



2/ = j J- ~ In 



5 ft /8a is fixed by the constraint of number conservation 
and is generally a complicated function of a and j3. The 
constancy of a\ D + ($\ D in one dimension suggests that 
we look for a similar condition in two dimensions. In the 
limit that the polarization is small relative to the filling, 
numerical results confirm that 



5_l 
8 a 



a 

1' 



(55) 



We note that /3 is nonzero for any finite a R - a single- 
particle Rashba interaction always polarizes the system. 
We may minimize Eq. (|46l) with a R strictly zero to obtain 
the critical J for symmetry breaking 



Jc = —t 

a a 



(57) 



which is consistent with Eq. (|4I|) in one dimension. Be- 
yond the approximation in Eq. (|55[) . we examine a correc- 
tion to Eq. ([57|) in Sec. V-B and find that it is small for 
moderate values of V. For a half-filled band in the unpo- 
larized state, electrons occupy all momenta with positive 
cos k x + cos k y , hence symmetry breaking is most likely 
to occur at half filling where a is maximized. We derive 
in App. C the following equation 

Jc K 

t 



(2tt) [ ki 



].(> 

768 



i- 8 
9216 



(58) 



0.8 



co. 
+ 



0.6 = 



0.75 



0.4 



0.2 



0.5 



0.25 



0.2 0.4 0.6 0.8 

Number of polarized electrons 



FIG. 4: a +/3 vs number of polarized electrons. Parameters 
are t = 1,J = 1,V = 1 and a R is varied to change the 
extent of polarization. Numbers 0.25, 0.5, 0.75 and 1 on the 
graph refer to the total number of electrons, a 2 + /3 2 can 
be approximated as a constant until the number of polarized 
electrons exceeds half the total number of electrons. 
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that is valid for small filling; J\(k F ) is a Bessel function 
of the first kind, of integral order 1. In this limit we 
see explicitly that J c is lowered when filling increases. 
For a R — and J < J c , we may express the bandwidth 
defined in Eq. (gSJ) and (J5TJ) as 



2t 



J 



(59) 



As J is increased, symmetry breaking prevents the 
bandwidth from reaching the critical value of aoV/2. 




FIG. 5: Analytic approximation of J/t vs a R /t. These are 
the values of J, a R and t consistent with a polarization of 0.1 
for fillings 0.25,0.5,0.75 and 1. 

We observe that a and /? are just the band filling and 
polarization weighted by different geometric factors in 
their respective integrals. For fixed a and /3, we may ask 
what are the range of J/t and a R /t that are consistent 
with this choice of polarization and filling. Eq. (|56|) gives 
us a simple linear relation between J/t and a R /t that we 
plot in Fig. (5) for various fillings. A Rashba-polarized 
system that is conventionally explained by a nonzero a R 
is indistinguishable from a system with a smaller (even 
zero) a R but larger J - there is a trade-off between a R 
and J which is characterized by the slope Typically, 
— /3 ~ 0(p) with p the number of polarized electrons 
normalized to a maximum of 1; this maximum occurs 
in the fully-polarized half-filled band. We explore the 
implications of Eq. (jS"6")) for Au(lll) in Sec. III-D. 



The analytic approximation in Eq. (|55[) is equivalent 
to the statement that the polarization is independent 
of V. We have shown that this approximation is only 
consistent if the resultant polarization is small relative 
to the band filling. We defer a discussion of the role of 
V beyond this approximation to Sec. V-B. 



D. Comparison with j4u(111) 

We investigate the range of values of t, V, J, a R 
that are compatible with the observed polarization of 
j4 it (111) surface states, which has a hexagonal Brillouin 
zone. The small filling suggests we may approximate the 
Brillouin zone as square. 



Au(lll) has a lattice constant of 2.89A. The Fermi 
momenta of the two Rashba bands are 0.153A" 1 
and O.mA^Hol]. The energy dispersion is fitted to 
E(k) = 1.82eV(/fc±0.0338) 2 . In the low-fc limit, the tight- 
binding dispersion reduces to E = — a(2 — k 2 /2) ± bk 
from which we infer a ps 3.64eV and b s» 0.123eV. 
a and f3 are evaluated by numerically evaluating 
Eq. (T4"3"l) and (|4"4")l assuming circular Fermi surfaces: 
a = 0.0707,/? =-0.00244. 

We define n (p) as the band filling (polarization) 
normalized to the range [0, 2]([0, 1]). Given the above 
values of a and (3, we numerically solve the self- 
consistent mean field equations to obtain n — 0.0363 
and p = 0.00525. This can be compared with the 
fractional filling and polarization obtained if we assume 
the experimental dispersion is embedded in a square 
lattice, i.e. n = tt(0.153 2 + 0.177 2 )/(2tt) 2 = 0.0364 and 
p = tt(0.177 2 - 0.153 2 )/(2tt) 2 = 0.00526. 

Since the polarization is small compared to the filling, 
we may employ the approximation of Eq. (|55p . From 
Eq. ([56]), J/t = 28.3 - 410 a R /t; this is plotted in Fig. 
(6). 




0.08 



FIG. 6: J/t vs a R /t in the model of Au(lll). 

The conventional understanding of ^4m(111) is equiv- 
alent to setting J = V = 0, hence a R /t — 0.069. We 
ask if the experimental data is consistent with a smaller 
value of a R . 
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(i) From Eq. (|5Tj) . the energy splitting is related to the 
Rashba constant and /3 by: a R + 0.00122 V « 0.123 eV. 
Reducing a H requires that V be unphysically large if the 
momentum-space polarization j3 is small. As per our 
discussion in Sec. III-B, the effect of V is maximized at 
half-filling; A«(lll)'s band is near-empty, hence V by 
itself is too weak to produce the required energy splitting. 

(ii) The energy-minimizing condition in Eq. (|56[) 
relates J, a R , t, a, (3: J/t — 28.3 — 410 a R /t. Reducing 
a R while keeping t fixed requires that J be unphysically 
large. This is because J most effectively polarizes a band 
at half- filling but the ^4m(111) band is nearly empty. 
One option to avoid enlarging J is to decrease both a R 
and t simultaneously. However, this option is costly in 
V as per our discussion in (i). 

We conclude that it is unlikely that J and V play a 
significant role in ^tu(lll). 



with the approximation of Eq. (|55[) . 



IV. CASE STUDY OF THE SPIN-SPLIT STATE 

In Ref. a parity- breaking spin-split Fermi-surface 
instability was studied with a reduced Hamiltonian that 
limited interactions to only nearest-neighbor exchange 
(J) of electrons with antiparallel spin. The Hamiltonian 
in Ref. [3] is equivalent to the particular choice V — J 
and At = in the more general reduced Hamiltonian 
considered in this paper; we demonstrate this in App. 
D. 



A. One Dimension 



E. Effects of At and J' 

The near-identical structures of Eq. (|2U|) and Eq. (|2"Tj) 
(except for minus signs) suggest the opposing roles played 
by J' and V. One can show that the only effect of J' is 
to cancel V. We may replace 



V 



K« = V - J' 



(60) 



in the preceding analysis. On polarization, J' helps to 
increase the bandwidth but reduces the energy splitting 
- its effects are self-cancelling. 

The effect of At on the energy spectrum is 

E n (k) — > E v (k) +2At(a + n(cosk x +cosk y )^j. (61) 

At promotes polarization in two ways: 

(i) The hopping parameter t is renormalized to 
t eft = t — At n. This reduces the kinetic energy cost of 
polarization. At is the only nearest-neighbor interaction 
that breaks particle-hole symmetry on a square lattice - 
this is reflected on the dependence of t 

eff on 71. 

(ii) On polarization, a decreases and the energies of all 
electrons are reduced by the same amount, independent 
of k. 



The effect of At on the ground state energy is 

^— '- — > ^-!- +2Atna. 
N N 

The energy minimizing values of a and ft are related by 



1 2 , A , 

a a v ' 



(62) 
dby 
(63) 



The spin-split state is equivalent to the Rashba-likc 
state in one dimension. In Fig. (7-a), we depict a 
spin-split state in which the spin is quantized along a 
direction perpendicular to the electron's momentum; 
in Fig. (7-b) we show the equivalent Rashba-polarized 
state. As a result of this equivalence, they share the 
same dispersion (Eq. ([36])') and ground state energy (Eq. 
1(57)) with a R = 0. 



(a) 



(b) 



FIG. 7: (a) The spin-split state in one dimension. Spin up 
and spin down electrons are shown divided by the momentum 
axis, (b) The Rashba state in one dimension. Opposite 
Rashba polarizations are shown divided by the momentum 
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B. Two Dimensions 

The equivalence between the spin-split and Rashba- 
like states in one-dimension does not generalize to two 
dimensions. 

We begin with the general Hamiltonian in Eq. 
Since the spin quantization direction is arbitrary in the 
spin-split state, we choose for convenience that it lies 
along z. We decouple the electron-electron interaction 
with the mean fields 

n k = rifct + n*4 = {cj^cvf) + (cj^Cfcj,) and (64) 
A fc = rifc t - n*4 = (c^Cfef) - (c^Cki) (65) 

The spin-split state satisfies the symmetry con- 
straints: rik = Afc = — A_fe. Apply- 
ing these symmetry constraints and the identity 
J2a { n kc n k > ±a ) = {n k n k > ± A fe A fe ')/2, we reduce this 



Hamiltonian and keep only H t ,Hj , if, 



(9=0) rr(9=fe'-fc) 



and Hf, =0) . 



We define geometric factors 



v 



= Jr X!( C0S kx + C ° S ^ nk ' 
k 

= 7j^2 { sink x + sinfej,) A fc 



(66) 
(67) 



and use the convention: a = +1( — 1) for spin up (down) 
The energy dispersion is 



E a (k) = - 2( t cff - i(2J - V cS )a ) (cosk x + cosk y , 



a \ g VcS 7 I ( sin kx + sin J 



(68) 



The ground state energy is 



^ = -2t cS a + ^Ja 2 -^V cS {a 2 + 1 2 ). (69) 



The mathematical structure is very similar to that 
of the Rashba instability (c.f. Eq. (gSJ) and flj5]) ). 
The difference lies in geometry: (i) a is different for 
the Rashba and spin-split states because the electrons 
occupy different regions in momentum space (ii) (5 and 
7 have different geometric weights in their respective 
integrals. 

We observe that the two-dimensional dispersion (Eq. 
(|68p) decouples into two one-dimensional dispersions as 



follows Q3]: 

Ea{k) =J2 ( - 2 (*eff- \(2J-V eS )a v )cosk v 



- a(2V e s'y v )smk l 
: 2]yX] cosfel/ ™ fc ' 

k 

: ^^sinfc„ Afc. 



(70) 
(71) 

(72) 



Here the index v runs over the {x, y} directions. If there 
is an energetic gain to polarizing in the x direction, 
the symmetry of the lattice permits the same gain for 
a polarization in the y direction. Hence, there is a 
spontaneously broken symmetry in which the two spin- 
differentiated Fermi surfaces are displaced in opposite 
directions along either x-y diagonals. One example 
of the Fermi energy contours of the spin-split state is 
illustrated in Fig. (1-a). 

This suggests we may think of the two-dimensional 
spin-split state as being composed of many one- 
dimensional systems aligned parallel to the x-y diagonal. 
For each one-dimensional system, number conservation 
ensures that the ratio a 2 + j 2 is a constant upon 
polarization; this is demonstrated in Eq. (|52p and (|53|) . 
With the identity, a — 4a„,7 = 47^ for any v, we 
have proven that the geometric combination a 2 + y 2 is 
rigorously a constant independent of polarization. An 
alternative proof is presented in App. E. 

Minimizing (H) while keeping a 2 + j 2 constant, we 
obtain for the spin-split state 



a 



2 tpff 



The critical J for symmetry breaking is 

_ 2 t c ff 

'J c — ■ 

The dispersion simplifies to 

toS Vrff 



(73) 



(74) 



J 



{cosk x + cosky) 



1 



V cS ^(smk x + sinfc y ). 



(75) 



Eq. (|64p to (|68[) and the energy-minimization condition 
(E t (k)) FSt =(E^k')) FSi (76) 

may be solved numerically for the unknown factor 7. Al- 
ternatively, we may numerically evaluate a for the case 
of zero polarization and extract 7 from 



a 2 = a 2 +1 2 = ^ +1 2 . 



(77) 
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In the limit of small filling, we use a result from App. C 
to obtain 

7 = I sKk* Ji(m) 2 -(t^) 2 - (78) 
V. THE ROLE OF V 



However, we must consider the factor in these 
integrals. For the same set of parameters, there are 
generally more electrons polarized in the Rashba state 
because electrons polarize in all directions, while elec- 
trons in the spin-split state polarize only in the x-y 
diagonal direction. We support this claim numerically 
in Fig. (8). 



In our previous analysis of the Rashba-like (Sec. Ill) 
and spin-split (Sec. IV) instabilities, we have focused 
on the role of nearest-neighbor exchange J in driving a 
symmetry breaking but have said little about the role 
of nearest-neighbor repulsion V . In this Section we ask 
how V (or effectively V c s — V — J') affects these two 
states. We find that V e s is pivotal in deciding which of 
the two ground states are energetically favored - this is 
discussed in Sec. V-A. In Sec. V-B, we comment on 
how V e ff affects the Rashba-like polarization beyond the 
analytic approximation of Eq. (|55|) . 



A. Competitition between the Rashba-like and 
spin-split instabilities 

The Rashba-like instability with a R = is energet- 
ically compared with the spin-split instability. Except 
for a small correction that is explained in Sec V-B, they 
have the same critical J c = 2t/a . For convenience, we 
call this the first critical J cl . Which state is preferred for 
J > J c i is decided largely by a competition between the 
geometric combinations a 2 + j3 2 and a 2 + 7 2 ; the param- 
eter V plays a larger role than parameters t e g and J. To 
understand this, consider the difference in energies of the 
Rashba and spin-split states: 



6(H) = S(H t , 



S(H V 



(79) 



The first term in this difference is ~ O(Sa) 2 because 
(-Ht eff , j) is minimized with respect to a for the spin-split 
state. The second term is ~ 0(Sa, 5/3) and is decisive in 
determining where the spin-split and Rashba states are 
degenerate. 
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FIG. 8: The absolute number of polarized electrons vs J/t for 
the Rashba (labelled R) and spin-split (labelled SS) instabili- 
ties at different band fillings: n = 0.25, 0.5, 0.75. Parameters: 
t — 1, Vcft = 3. For all fillings, the critical J c \ for symme- 
try breaking from the paramagnetic phase is nearly equal for 
both Rashba and spin-split states. 



Small filling 



We first examine the geometric difference between (3 
(Eq. ([HI) ) and 7 (Eq. (|o7j)). As a first approximation, 
this difference suffices to explain qualitatively the behav- 
ior of a 2 + f3 2 relative to a 2 + 7 2 . For any momentum, 
the geometric weight in the integral of j3 is always less 
than or equal to the weight in 7: 



sin 



' k y < sin k x + sin k y 



(80) 



In one dimension these weights are trivially equal, hence 
the spin-split and Rashba states are exactly equivalent 
(see Sec. IV-A). In two dimensions, Eq. (|8T)|) naively 
suggests that /3 < 7. 



For small filling, the fractional difference in the 
number of polarized electrons is small - the effect of 
is outweighed by the difference in geometric weights. 
Hence /3 < 7. This explains why a 2 + f3 2 is a mono- 
tonically decreasing function of J/t for small fillings; 
we show this in Fig. (9-a). Since a 2 + f? < a 2 + j 2 , a 
comparison of the energies in Eq. (pfrJ|) and (IMl) implies 
that the spin-split state is favored for small filling if 
V c ff = V - J' > 0. We demonstrate in Sec. VII that 
a large positive V e s is necessary for either spin-split or 
Rashba state to exist at all, hence we predict that the 
spin-split state exists only for small filling. 
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0.167r 



a 2 +y 2 



rection ~ O(5o) 2 due to 5(H toii .j). We recall from Eq. 
P3"| that a, the relevant geometric factor in (Ht eBt j), 
is the negative sum of the kinetic energy; a is smaller 
for the spin-split state because it polarizes in only one 
direction and incurs a larger kinetic energy cost. Hence, 
{Ht ea ,j)ss < {H teB ,j)n and the effect of t oS and J 
through 6{Ht etli j) is to reduce J C 2- 



0.164 
4.8 



5.4 
J/t 

(a) 



We may ask how the magnitude of V e g affects J C 2- 
In the limit of large V c g, 5(Hy eS ) ^> 5(Ht etlt j) hence 
Jc2 Ji- The limit of V c r — ¥ is unphysical for two 
reasons: (i) we note from Eq. (ITS! that the spin-split 
state has a bandwidth and energy splitting proportional 
to V e g (ii) in Sec. VII we show that both spin-split and 
Rashba states are unstable against ferromagnetism for 
small Veg. In Fig. (10), we plot J c i,J C 2 and Jj as a 
function of band filling for V e g = 3. Above the critical 
filling and for J cl < J < J c2 « Jj, a 2 + /3 2 > a 2 + 7 2 so 
the Rashba phase is favored; for J > J c2 , the spin-split 
phase is favored. 




2. Phase transition above a critical filling 

Above a certain critical filling (n c — 0.45), an interest- 
ing geometric feature of a 2 + (3 2 leads to the possibility 
of a transition between Rashba and spin-split states at 
a second critical J C 2, where these states are degenerate 
for a second time. In Fig. (9-b), we observe that for 
quarter-filling, a 2 + j3 2 initially rises with increasing J, 
then eventually falls and intersects a 2 + 7 2 when J = Jj. 

As we have argued through Eq. ([79]) . the spin-split 
and Rashba states are degenerate near this intersection 
Jj. J C 2 would lie exactly on Ji if not for a small cor- 



FIG. 10: Phase diagram for the Rashba (R), spin-split (SS) 
and paramagnetic (P) states as a function of band filling and 
J/t. V e g = 3. a labels J c i; b labels J C 2\ c labels Ji, the 
intersection of the geometric factors a 2 + and a 2 + f3 2 . 

We understand the behavior of a 2 + (3 2 and a 2 + -f 2 
in the following manner: 

(i) In the limit of large J/t, the polarizations of both 
Rashba and spin-split states saturates to a maximum 
(p — > n) and the fractional difference in the number of 
polarized electrons diminishes. The geometric weights 
dominate the integrals in j3 and 7, hence /3 < 7. 



14 



(ii) In the limit of small (J — J c i)/t, the Rashba 
state has significantly more polarized electrons than the 
spin-split state, as we may observe in Fig. (8). This 
allows Afc to dominate the integrals in j3 and 7. The 
effect of Afc becomes more pronounced as n — > 1 because 
the fractional difference in the number of polarized 
electrons increases with the band filling. Consequently, 
there exists a critical filling n c = 0.45 above which j3 > 7. 

From (i) and (ii), we deduce that there must be an 
intersection J[ > J c i where a 2 + (3 2 = a 2 + j 2 . In addi- 
tion, (ii) implies that as the band filling is increased, the 
range of J/t for which /3 > 7 also increases, hence both 
Jj and J C 2 increase. This increase is quite dramatic, as 
we show in Fig. (10). The spin-split phase is effectively 
squeezed out of existence by the Rashba phase for large 
filling because J C 2 ~> t. 

The introduction of the single-particle Rashba in- 
teraction [a R > 0) changes the phase diagram in Fig. 
(10) in two ways: (i) the spin-split phase will shrink in 
parameter space (ii) there is no longer a spontaneous 
broken symmetry for the Rashba state; the paramagnetic 
phase becomes Rashba-polarized. 



breaking 



V cS fS(a 2 



5 a 



(82) 



We have denoted the critical J of the spin-split state 
as J^ s — 2t c s/a - The second term is a correction 
to the analytic approximation of Sec. III-C. From 
our discussion in Sec. V-A, we deduce that J c > J^ s 
(J c < J^ s ) for band fillings less than (more than) 0.45. 
The ratio of (8(a 2 + f3 2 )/5aj o is fixed by the constraint 
of number conservation and is typically much less than 
a Q for Rashba-like systems - J c w J* s for moderate 
values of Kff- 



VI. GENERAL CLASS J GROUND STATES 

In this Section we present the formalism for the 
description of a general Class J ground state. We 
consider the energetic stability of Class J states other 
than the spin-split and Rashba-like states and comment 
on various symmetries such as parity and time-reversal. 



The Role of V e a in Rashba-like states 



Formalism 



In our previous discussion of the Rashba-like state in 
Sec. III-C, we have made an analytic approximation 
in Eq. (l55j) that is equivalent to claiming that the 
polarization is independent of V c g. We have shown that 
this approximation is only consistent if the resultant 
polarization is small relative to the band filling. We now 
go beyond this approximation and ask how V e g changes 
the Rashba polarization. 

We consider an infinitesimal variation to the ground 
state energy in Eq. (|4"o| due to an increase in polarization 



N 



{ J a — 2 t c ff)6a 



VcS xf 2 
—6{a 



/? 2 



(81) 



VeH assists polarization if 6 (a 2 + (3 2 ) is positive. From 
our discussion in Sec. V-A, we know that the sign of 
(5 (a 2 + /3 2 ) is always negative for band fillings below 
the critical level of 0.45 - we conclude that V e g discour- 
ages Rashba-like polarization for bands with small filling. 



For fillings above the critical level, we note from Fig. 
(9-b) that the slope of a 2 + f3 2 levels off at some critical 
polarization which we call p c . V e & assists (suppresses) 
polarization for p < p c (p > p c ); in the limit of large V e g, 
the polarization will approach this critical value. 

From Eq. (|81l) we derive the critical J for symmetry 



A general Class J ground state may be character- 
ized by their spin structure ip(k) (defined in Eq. ((SJ)) 
and polarization structure A(fe) (defined in Eq. (fTU)) ). 
The various possible ground states are distinguished by 
the component of the interacting Hamiltonian that de- 
scribes nearest-neighbor repulsion; for example, see Eq. 
(|24l) . From Eq. ([24]) we may separate the amplitude for 
nearest-neighbor repulsion 

y-^ v e i(k-k'ys 

s 

= 2 V I sin k v sin k' v + cos k v cos k' v J . (83) 

We consider only the component of Eq. (|24p that is de- 
pendent on ip and A 



N 



TEiflv-E^ sinfc. 



l^A fcX (fe)sinfc^ j. (84) 



In arriving at Eq. (|84|) we have dropped the sec- 
ond term proportional to cosfc„ in Eq. ([55]) ; this is 
necessary to satisfy the symmetry constraint in Eq. (|12l) . 
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We define geometric factors 

C, =— A ( k ) C(*0 sin k v and (85) 

k 

Xu=^J2Mk) X {k)smK. (86) 

k 

Eq. (|55|) and (|56")l suggest we may interpret („ an d 
Xjz as generalized Fourier coefficients of the harmonics 
£(fc)sinfc„ and \(k) sink^ respectively. The change in 
the energy dispersion due to £ and \ is 

4 C ' X) (k) =-nV^2 K U{k) ( v + X (k) X u\ ■ 

(87) 

Including the effects of J', the ground state energy be- 
comes 

= ~2t eS a + \jc?- \ VeS (a 2 + 2 £ (g + . 

(88) 

B. Energetic Considerations of other Class J 
Ground States 

In Sec. Ill, IV and V we chose to analyze the spin-split 
and the Rashba-like states because they are the simplest 
examples of Class J ground states. In this Section, we 
demonstrate that they are also the most energetically 
stable phases in their class. The general principle behind 
these claims is that anisotropy is energetically costly. 

In Sec. II we learned that the spin polarization 
as a function of momentum is given by ip{k) — 
£(fc) + ix(fc) = ( aX ) ~ i ( aV )- For simplicity of notation 
in this Section, we will omit the implicit dependence on 
momentum for all functions. Given a fixed spin struc- 
ture ip, Eq. (|87p informs us that A may spontaneously 
develop a symmetry in any of the four possible har- 
monics - {Csin/cj,, xsinfci,} - that satisfy the symmetry 
constraint of Eq. (fTTj) . It is also possible to combine 
two harmonics to obtain a more favorable ground state. 
A case in point is the spin-split state. Taking ip = 1, 
we have two possible harmonics - sin k x and sin k y 
- which correspond to a displacement of the Fermi 
surface in the x and y direction respectively. These 
harmonics are complementary in the sense that we may 
displace the Fermi surface along the x — y diagonal for 
a greater reduction in energy. Hence, the spin-split 
state has a polarization with the combined symmetry 
A ~ sin k x + sin k y . 

The spin-split example suggests that a state with the 
maximum number of complementary harmonics is likely 
to be energetically favored. We consider a general state 
that develops a symmetry in all four harmonics; we ask 



what conditions are necessary to keep all four. We define 
the reflection operations R x by [k x — > ~k x , k y — > k y ) 
and R y by (k y — > —k v , k x — > k x ). For an arbitrary 
band filling, the condition of no net magnetization 
(Eq. (fTTj) ) enforces certain conditions on how £ and x 
transform under R x and R y . The polarized spin ip A has 
six distinct harmonics {(^sin/cj,, x 2 sinfc„, Cxsin/c^} 
which must all sum to zero to satisfy Eq. (fTTj) . i.e. for all 
momenta tjj A — ¥ —ip A under inversion - there are only 
two possibilities, which we categorize into subclasses of 
Class J. 

The first subclass of ground states is defined to satisfy: 
(i) ( is even under both R x and R y (ii) \ is even under 
both R x and R y . As a result, A(fc) = — A(— k) - the 
Fermi surfaces are displaced relative to each other in 
opposite directions. The choice of 9 independent of 
momentum corresponds to the spin-split state. 

The second subclass is defined to satisfy: (i) C is even 
under R y but odd under R x (ii) \ is even under R x but 
odd under R y . Alternatively, the conditions (i) and (ii) 
describe an equivalent state if ( and x interchanged. 
The combined effect of (i) and (ii) is that the spin 
polarization rotates 2ir radians around the Fermi con- 
tour - it has winding number one. The polarization A 
has (i) two 's-like' harmonics {xsink Xl (^smk y } which 
are even under R x and R y (ii) two 'd-like' harmonics 
{(^sinfcy, x sin fez } which are odd under R x and R y . 
Unlike the first subclass, the second subclass has two 
competing symmetries and will favor one and not the 
other. We note that both 's-like' and 'd-like' instabilities 
satisfy ip A — > —ip A under inversion; we elaborate on 
this in the next Section. 

The simplest example of the second subclass is 
i-' = /fc/|/fe|i fk — sinfca; ± zsinfcj,; this has three distinct 
harmonics {sm 2 k x /\f k \, sin 2 fc y /|/ fc |, sin k x sin k y /\f k \}. 
The first two correspond to a relative expansion of 
one Fermi surface in the x and y directions; they are 
complementary in the sense that the Fermi surface may 
expand in both directions with a combined symmetry 
A ~ y/ sin 2 k x + sin 2 k y - we call this the Rashba-like 
phase. The third harmonic corresponds to an expansion 
along one x — y diagonal and a contraction along the 
other; this is not complementary with the symmetry 
Y^sin 2 ^ + sin 2 fc !/ . 

This general discussion suggests a great many pos- 
sibilities among the two subclasses; however, we now 
argue from energetic considerations that only two really 
matter. The Rashba-like and the spin-split phases 
are special because they have the minimum amount 
of anisotropy needed to reduce the magnetic moment 
of a polarized state to zero. The Rashba-like phase 
has the maximum isotropy in charge space and the 
next-best isotropy in spin space and vice versa for the 
spin-split phase. In the continuum analogy, the Rashba 
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phase corresponds to an s-wave charge polarization and 
a p-wave spin polarization, while the spin-split phase 
corresponds to a p-wave charge polarization and an 
s-wave spin polarization. 

A highly anisotropic harmonic {Csinfc„, ^sinfcj,} 
encourages an anisotropic polarization that is penalized 
by high kinetic energy (i.e. greatly reduced a 2 ). In 
addition, the 'surface tension' due to H t ensures that the 
overlap between the polarization A& and the harmonic 
is weak - the generalized Fourier coefficients {£ u , \v} 
are diminished. In combination, a 2 + 2J2i> (Cv + xt) 
for these anisotropic ground states are generally smaller 
than that for the spin-split and Rashba-like phases; 
from Eq. (|88p. we conclude that the Rashba-like and 
spin-split states are the most favored. 

We may ask physically why the Rashba-like and the 
spin-split states are so close in energies. In the spin-split 
phase, electrons with neighboring momenta have spins in 
perfect parallel alignment, hence this state benefits from 
the most reduction in energy per polarized electron 
due to the exchange interaction V e e; this exchange in- 
teraction is explained in Sec. II. However, the spin-split 
state polarizes in one direction and consequently has 
fewer polarized electrons than the Rashba-like phase. 
The Rashba phase has a spin structure that rotates 
2tt as we go around the Fermi contour; the spins of 
electrons with neighboring momenta are not perfectly 
parallel to one another, hence the energy reduction per 
polarized electron is not as large. The net effect is that 
the Rashba-like and spin-split states are almost equally 
favored by the exchange interaction V e e. 



2. Condition for Time-Reversal Invariance 

For time-reversal invariance to hold, we require that 
for all momenta tp A — > —yt> A under parity inversion, 
i.e. any polarized electron always has a time-reversed 
partner with opposite momentum and spin polarization. 
In Sec. VI-B we showed that this condition is satisfied 
for the first subclass and both 's-like' and 'd-like' in- 
stabilities in the second subclass; they share the crucial 
condition that £ x is even under inversion. 

If we relax this condition, it is still possible to obtain 
a state with no net magnetization if £ = £(k x ) is even 
under R x and x = x(kx) is odd under R x . A has two 
competing harmonics: (i) Q{k x ) sin k y is 'p-like' (ii) 
x(k x ) sinfcj, is 'd-like'. For the favored 'p-like' instability, 
the y-th component (— %A) of the polarized spin has 
the symmetry sin (2 k x ) sin k y which is invariant under 
inversion, hence time-reversal symmetry is broken. 
This general discussion suggests that time-reversal 
breaking states in Class J have at most one har- 
monic and are not stable relative to multiple-harmonic 
states like the spin-split and Rashba-like states. The 
simplest state that breaks time-reversal symmetry is 
^ = e- jfe *, Cx = sin(2fc x )/2. 



VII. COMPARISON WITH THE 
FERROMAGNETIC INSTABILITY 



C. Comments on Various Symmetries 

In this Section we ask if Class J grounds states (i) 
break parity (ii) preserve time-reversal symmetry. We 
provide a proof for the first and a necessary condition 
for the second. 



The ferromagnetic instability has a net magnetic 
moment and is not a Class J ground state. The ferro- 
magnetic ground state has the most isotropic 'polarized 
spin' - it corresponds in the continuum analogy to an 
s-wave charge and s-wave spin polarization. In this 
Section, we compare the ferromagnetic with the Class 
J instabilities. We single out nearest-neighbor repulsion 
V as the crucial interaction that allows Class J ground 
states to compete with ferromagnetism. 



1. Proof of Parity-Breaking 

At least one component of the polarized spin must 
have a harmonic of the form £ 2 sin k v or x 2 sin k v . On 
inversion, harmonics of this form changes sign. This 
implies that the spin of the polarized electrons are not 
invariant under inversion, hence Class J ground states 
break parity. 



We begin with the Hamiltonian in Eq. ([J} and 
decouple the electron-electron interactions with the 
mean fields defined in Eq. (|64|) and (|65|) . In contrast 
with the spin-split state, the ferromagnetic state satisfies 
the symmetry constraints: = n-k, A& = A_^. 
This difference changes the mean-field Hamiltonian 
significantly. Following a previous work by one of us 
[H, we keep terms: H t , H{j =0) , Hj, H^ =h '- k) , H£= 0) 
a,nd H ( j q r 0} . 
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We define the geometric factors 

a =— ^2(cosk x + cos k y ) n k , (89) 

k 

£ =— y^(cos k x + cos ky) Afc and (90) 
fc 



(91) 



The energy dispersion is 

E ka =( -2t + 2Atn+ -(2J + J'-V)a\ ^ cos k u 

\ ~ / y 

-cr ( - ( J' + V)£ + 2 At k J ^2 cos fc ^ 

-a(2At£+ ^(4J + U)k\ +2Ata (92) 
and the ground state energy is 
^ =-2(i- Atn)a+i(2J + J' - U)a 2 

- ^(4J + U)k 2 - + V)£, 2 -2A*k£. (93) 

Unlike the Rashba and the spin-split states, the fer- 
romagnetic state is characterized by a spin-dependent 
bandwidth and energy baseline. 

There are more channels through which the ferromag- 
netic state may gain an energetic advantage. It is no 
surprise that within the tight-binding Hubbard model, 
the ferromagnetic instability is favored for a large range 
of parameters. 



A. One Dimension 

In Sec. IV, we showed that the Rashba and spin-split 
states are equivalent in one dimension - for convenience, 
we call this the RSS state. In this Section, we ask: for 
what range of parameters is the RSS state favored over 
the ferromagnetic state in one dimension? 

We parametrize the ferromagnetic polarization with di- 
mensionless quantity e, which is defined: n ka = 1 f° r 
k € {— k F — ere, +k F + ere}. Minimizing the ground state 
energy with respect to e, we obtain an implicit equation 
for the ground state polarization e* : 



2tt t sin k F — 4 At k F sin k w — 4 At cos k w ] sin e* 



U + 2Jj e* + l^At cos k F j e* cos e* 
- I 2Jsin 2 fc F + V cosk F + J' ) sine* 



The Hmit of vanishing e* determines the critical J c for 
the onset of ferromagnetism ,25] 

J c = ^7r t sin k F — 2 At(k F sin k F + 2 cos k F ) 

-j-Vcos2k F -J'j/(l + 2sm 2 k F y i (95) 



from Eq. (|37| and (|60|) we extract the critical J c for the 
RSS state 



Jc 



7T t - 2 At k F 
2 sin k F 



(96) 



(94) 



If J is greater than both J c and J c , ferromagnetism typ- 
ically wins because it gains an extra energetic advantage 
through At,U and J' (c.f. Eq. gBJ),®,®). 

We analyze the effect of the various Hubbard parame- 
ters: 

(i) In Eq. (|9"5")l . we note that V is multiplied by cos 2k F 
- a large V increases J c in the range k F £ [tt/4, 3n/A]. 
This allows the RSS state to compete with the fer- 
romagnetic state. This is illustrated in Fig. (11) for 
parameters t = 1, At = 0, U = 0, J' = 0, a R = 0. The 
RSS state exists in the region where J c > J > J c , which 
is most likely to occur at half-filling and extends to 
quarter-filling in the limit of large V. For small V , J c 
never intersects J c (compare curves a and b in Fig. (11)) 
and ferromagnetism is always favored. 

(ii) Positive parameters U and J' tend to favor 
ferromagnetism by reducing J c . This narrows the region 
in which the RSS phase is favored. 

(iii) The presence of At breaks the particle-hole sym- 
metry of the phase diagram. It improves (reduces) the 
stability of the ferromagnetic phase at small (near-full) 
filling respectively. 

The introduction of a nonzero a R will Rashba-polarize 
the paramagnetic phase and shrink the spin-split phase; 
if a R j3 is comparable to the energetic gain for ferro- 
magnetism, the ferromagnetic phase will diminish as well. 



B. Two Dimensions 

In one dimension, we emphasized the role of nearest- 
neighbor repulsion V in allowing the Rashba-like and 
spin-split states to compete with ferromagnetism near 
half-filling. In this Section we ask if this conclusion 
generalizes to two dimensions. 

We proved in Sec. III-B that the energy-minimizing 
condition is equivalent to equating the single-particle dis- 
persions (Eq. (1921) ) of spin up and down electrons at their 
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FIG. 11: Phase diagram for the RSS, ferromagnetic (FM) and 
paramagnetic (P) states as a function of the Fermi momentum 
and J/t; one dimension; parameters: t = 1, At = 0, U = 
0, J' = 0. a labels the critical J c for ferromagnetism with 
V = 5; b labels J c for ferromagnetism with V = 1; c labels 
the critical J c for the RSS state. 



respective Fermi energies. We define these Fermi energies 
as e% and £^ and perform said operation 



4Ai £ 



4 J+ r ) /, 
1 



At_ 2J + J'~V 
— n a 

t 4t 



J' + V At 



At 



t 



(97) 



As is well-known in a mean- field treatment, ferromag- 
netism is greatly enhanced where the density of states 
diverges. In a one-dimensional nearest-neighbor tight- 
binding model, this divergence occurs for near-empty 
or near-full bands; in two dimensions the divergence 
occurs at half-filling. In contrast to one dimension, the 
two-dimensional Rashba-likc and spin-split states cannot 
compete against ferromagnetism at half-filling. 

We expect that the Rashba-like and spin-split states 
are stable if a large nearest-neighbor repulsion V sup- 
presses ferromagnetism; we know from Eq. (|98[) that 
this only occurs when the condition a > e p g(e F )/t is 
satisfied. In one dimension, we showed in Sec. VII-A 
that this condition restricts the band filling to lie in the 
range [1/4, 3/4]. In two dimensions, the density of states 
is approximately flat away from the divergence at half- 
filling. We approximate g(e) — 1 / (8t) and simplify Eq. 
to obtain 

/ J' 2 

J c = ( t - 4U - 2 At (2 - n) - — (l + (l - n) ) 



V 
~2 



l-3(l-nT 



2- (1-n) 



(100) 



This result was first derived by one of us [25|. Within 
this approximation, we find that V destabilizes ferro- 
magnetism for band fillings in the range [0.423,1.577]; 
this range is slightly larger than that in one dimension. 
At half-filling, the approximation breaks down and 
ferromagnetism dominates. We have learned in Sec. 
V-A that the Rashba-like phase likely outcompetes the 
spin-split state for band fillings larger than 0.45 or less 
than 1.55; this phase diagram is plotted in Fig. (TTU1) . 
As we increase the band filling from 0.423, we predict a 
transition from the spin-split state to Rashba-like state 
and then to ferromagnetism near half-filling. A more 
general treatment than this Section would also allow for 
antiferromagnetism at half-filling. 



We define the density of states g{e). In the limit of zero 
polarization, we find an expression for the critical J 



±g{e F ) + f t )J c 



1 - U g(e F ) + ^ ^4 e F g(e F ) - n 
V ( slg(e F ) \ J> 

+ u( a °-—r- --t {0 ' 



t 



(98) 



For a square lattice with only nearest-neighbor hopping, 
there is a Van Hove singularity at half-filling. Since the 
singularity in g(e) is logarithmic, the dominant terms in 
in Eq. (96) the limit e F — » are 



(99) 



g(e F )[ J c + U)=l 



VIII. CONNECTION WITH FERMI LIQUID a 
AND f3 PHASES 

Within a tight-binding Hubbard model with on- 
site and nearest-neighbor interactions, one of us first 
proposed that electron-electron interactions may sponta- 
neously produce a spin-split phase [lij]. In this paper, we 
demonstrated that the Rashba phase may be produced 
through a similar mechanism. Wu and coworkers have 
also explored the spontaneous generation of the a and 
/3 phases in a Landau- Fermi liquid framework [15]. The 
names of these phases should not be confused with the 
geometric factors defined in Eq. (|43|) and (|44)l . The 
spin and momentum-space polarizations of the spin-split 
and Rashba-like phases are equivalent, in the continuum 
limit, to that of the a and j3 respectively. In this Section, 
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we compare the formalisms and the conclusions derived 
from both models. 

In the Landau Fermi liquid framework, Wu et. al. 
showed that the a and (3 phases result from a Pomer- 
anchuk instability |l| in the I — 1 angular momentum 
channel and is driven by the phenomenological Landau 
parameter i* 1 " . The interaction part of their phenomeno- 
logical Hamiltonian may be simplified to 



(<z) 



with the spin current operator defined as 

Q AIb (9) = E c U, Q ^^. 



(101) 



(102) 



The coupling fi(q) is assumed to be a Lorentzian cen- 
tered around q = 0. They defined /" (0) = f£ and de- 
coupled this interaction with a q — mean field 



fib 



N 



k 



(103) 



We define the density of states at the Fermi surface 
as N(0). The Pomeranchuk instability occurs at 
JV(0)/f < -2. 



The Rashba ground state is achieved by choosing n 



fib 



Aib 



The resulting mean-field Hamiltonian is diagonal- 



ized by the Rashba basis: 



ffint^sgn^)^— f 



hij 



(104) 



This results in an effective spin-orbit field that Zeeman- 
splits the Rashba energy bands. We compare this 
phenomenological Hamiltonian with Eq. (j45|) and (|60l) . 
To cross over from the tight-binding to the Fermi liquid 
formalism, we limit interactions in the tight-binding 
model to a thin shell around the Fermi surface and take 
the continuum limit. By comparing the energy splitting 
of the Rashba bands in both formalisms, we may identify 
the phenomenological coupling /f with — V e gk F in two 
dimensions. 

In the tight-binding Hubbard model, the effective 
spin-orbit field is due to an exchange interaction in 
momentum space due to off-site repulsion V; there is 
also a cancelling effect by the pair-hopping J'. However, 
these same processes also increase the bandwith in a 
manner that suppresses polarization. In the continuum 
limit, we showed in Sec. V-B that the net effect of V e g 
is to suppress Rashba polarization. 

We reach a similar conclusion in the comparison of the 
a phase with the spin-split state. The phenomenological 
model postulates an effective Zeeman field that is 



necessarily accompanied by an increase in bandwidth 
within the tight-binding Hubbard approximation. We 
have shown in Sec. IV that these two effects exactly 
cancel and that V e f{ neither promotes nor suppresses the 
spin-split instability. We conclude that the validity of the 
phenomenological model must lie outside the physical 
approximations that are assumed in the tight-binding 
Hubbard model. Note that within the Fermi liquid 
analysis the a phase in two dimensions is predicted to be 
favored over the /? phase when the second derivative of 
the density of states at the Fermi energy is positive[l5j], 
which is the case in the band structure considered here 
for all band fillings (except half-filled) . Instead, in the 
lattice model we found that for band fillings above 0.45 
the Rashba phase is favored (Fig. 10). 

We have found that for Class J ground states, J is 
the greatest driving force for symmetry breaking. It is 
instructive to translate H { J =0) to the language of Pomer- 
anchuk instabilities [l|. In a basis J7, we decompose an 
arbitrarily-shaped Fermi surface into its Fourier compo- 
nents 



k n F {0) = K 



Ski 
V2 



.(U 



(105) 



We define k° as the Fermi momentum before symmetry 
breaking. From Eq. (O, the change in kinetic energy due 
to a deformation of the Fermi surfaces is 



N 



1 

47T 



l,Q 



(Ski 



(106) 



We compare this to the change in the energy of Hj 
From Eq. (|2"3"|). 



(9 = 0) 



8{Hf 



N 



(107) 



As we argued in Sec. II, a simplified Hamiltonian with 
only parameters t and J exhibits symmetry breaking in 
an arbitrary basis ft. The critical J is 4-7ri/(fc°) 2 which 
agrees with Eq. <(57j) and ([74]) to 0(k° F ) 2 . In addition to 
the arbitrariness of the basis, Eq. (I107j) also suggests 
that all Fourier modes I are unstable. We gain specificity 
when we take into account the effect of the lattice and 
of V e s; as we demonstrated in Sec. VI, the ground states 
with the least anisotropic 'polarized spin'- the spin-split 
and the Rashba-like phases - are the most favored. 



IX. DISCUSSION 

In this paper we have singled out a class of symmetry- 
breaking ground states with three properties: (i) 
electron crystal momentum is a good quantum number 
(ii) these states have no net magnetic moment (iii) 
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their distribution of 'polarized spin' in momentum space 
breaks the lattice symmetry. Examples of this class 
of ground states include the spin-split state and the 
Dresselhaus/Rashba spin-orbit coupled states. Since 
the Rashba, Dresselhaus and helicity spin-orbit coupled 
states are degenerate under Coulomb interactions, 
we just call them Rashba-like. Employing a mean- 
field approximation, we analyzed these ground states 
within a tight-binding single-band Hubbard model 
with all on-site and nearest-neighbor matrix elements 
of the Coulomb interaction. The relevant matrix ele- 
ments were found to be nearest-neighbor exchange J, 
pair hopping J' and nearest-neighbor repulsion V\ At 
plays the minor role of reducing the hopping parameter t. 

Because of these symmetry constraints, these ground 
states lower their energy most effectively through the 
nearest-neighbor exchange interaction J, hence the 
name Class J. Though J has been proposed to lower 
the exchange energy of localized electrons with parallel 
spin [23], the effect of J in momentum space is through 
the direct energy (26[. We have found that J favors 
a separation of electrons in momentum space, i.e. 
polarization. This conclusion is independent of the 
choice of basis, hence J alone does not energetically 
distinguish between different ground states. In the 
broken-symmetry phase, a mean-field decoupling of 
Hj reveals that J lowers the energy by expanding the 
bandwidth of the single-particle dispersion. 

The other relevant interactions are nearest-neighbor 
repulsion V and pair- hopping J'; they act in the 
combination V e s = V — J' to influence both the direct 
and exchange energies. In the broken-symmetry phase, 
V e ff (i) increases the direct energy by compressing the 
single-particle bandwidth (ii) reduces the exchange 
energy by splitting the dispersion of electrons with 
opposite polarizations. We have found that for ground 
states with highly anisotropic 'polarized spin', the 
advantage gained in exchange energy is outweighed by 
the cost in direct energy due to V c g - these states are 
energetically unfavorable. Among Class J ground states, 
the most favored are found to be the spin-split and 
Rashba-like states - they possess the minimum amount 
of anisotropy needed to eliminate their net magnetic 
moment. In the continuum analogy, the Rashba-like 
phase corresponds to an s-wave charge polarization and 
a p-wave spin polarization, while the spin-split phase 
corresponds to a p-wave charge polarization and an 
s-wave spin polarization. 

We provided a physical explanation for why the spin- 
split and Rashba-like states are so close in energy. In the 
spin-split phase, electrons with neighboring momenta 
have spins in perfect parallel alignment, hence this state 
benefits from the most reduction in exchange energy 
per polarized electron due to V c g. However, the spin- 
split state polarizes in one direction and consequently 



has fewer polarized electrons than the Rashba-like phase. 
The Rashba phase has a spin structure that rotates 2-7T 
as we go around the Fermi contour; the spins of electrons 
with neighboring momenta are not perfectly parallel to 
one another, hence the exchange energy reduction per 
polarized electron is not as large. The net effect is that 
the Rashba-like and spin-split states are almost equally 
favored by V e s; they are only distinguished by band 
filling. 

In the case of a square lattice, we have found that 
the spin-split state is favored for small band fillings; 
above the critical filling of 0.45, we predict a symmetry 
breaking from the paramagnetic to the Rashba-like 
phase at J c \ and a second phase transition from the 
Rashba-like to the spin-split state at J C 2 > J c i- Since 
all nearest-neighbor interactions on a square lattice 
except A t are particle-hole symmetric, V Q g produces an 
effective Zeeman field that is extremized at half-filling. 
This suggests that differences between the spin-split and 
Rashba-like states are magnified near half-filling - in 
particular, the Rashba-like state has significantly more 
polarized electrons than the spin-split state, hence the 
exchange energy favors Rashba. For near empty (or full) 
filling, the difference in polarization is not as pronounced 
and as a consequence the spin-split state is favored. 

We have found that a Rashba-polarized system that is 
conventionally explained by a relativistic single-particle 
interaction (a R ) with an electric field is indistinguishable 
from a system with a smaller (even zero) a R but larger 
J (and V e g for larger band fillings). The effects of J 
and V e ff are found to be strongest at half-filling. We 
have compared our model with the Rashba-polarized 
surface states in Au(lll). Due to its near-empty band, 
we concluded that electron-electron interactions are too 
weak to explain the large energy splitting in Au(lll). 

We have compared the energetic stability of Class J 
ground states with ferromagnetism, an instability with 
the most isotropic 'polarized spin'; in the continuum 
analogy, ferromagnetism corresponds to an s-wave 
charge and s-wave spin polarization. Because it has a 
less restrictive symmetry constraint, we found many 
more channels through which the ferromagnetic state 
may gain an energetic advantage. The ferromagnetic 
instability is favored for a large range of parameters, 
as is confirmed by the abundance of ferromagnetic 
materials in nature. The crucial interaction that allows 
Class J grounds states to compete with ferromagnetism 
is nearest-neighbor repulsion V. The particle-hole 
symmetric repulsion is found to maximally suppress fer- 
romagnetism at half-filling; in contrast, we have shown 
that its effects on the direct and exchange energies of 
Class J states are self-cancelling. Hence we expect 
that the Rashba-like and spin-split states are stable in 
systems with large V and only for a certain range of 
band fillings. From a mean-field treatment on a square 
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lattice, this range is found to be between quarter- to 
three-quarter-filling in one dimension; in two dimensions 
we argue that this range is enlarged to [0.423,1.577] 
with an exception at half-filling, where there is a Van 
Hove singularity. 

In this paper our analysis was done on a square lattice 
assuming only nearest-neighbor hopping terms in the 
single-particle Hamiltonian. As a result, all electron- 
electron interactions with the exception of correlated 
hopping At are particle-hole symmetric - it is not 
surprising that the effects of J and V e ff are maximized 
at half-filling. The conclusions that are specific to band 
filling will change (i) if next-nearest-neighbor hopping 
terms are relevant or (ii) if we analyze a different lattice 
(e.g. honeycomb, triangular). In particular, we expect 
that the regions of stability for the different Class J 
and ferromagnetic states will shift. For example, if we 
include next-nearest-neighbor hoppings on the square 
lattice, the Van Hove singularity moves away from 
half-filling - this may stabilize the Rashba-like and 
spin-split states at half-filling. In the three-dimensional 
generalizations of the spin-split and Rashba-like states, 
we believe our conclusions with regards to band filling 
are applicable to simple cubic and BCC lattices. 



X. APPENDIX 

A. Rashba tight-binding approximation 

We formally derive the Rashba tight-binding Hamil- 
tonian discussed in Ref. [28(. The free-space Rashba 
spin-orbit term is 

Hr = a R (a x k y - a y k x ) (108) 

for a single electron. We generalize this to a lattice as 
follows 

H R = ( Yk y C ko (^) , °k<y' ~ X k, cL (°y) , c ka' ) ■ 
koo' v 7 

(109) 

Here we have neglected band indices. 

We Fourier transform Eq. (|109|) and employ the sym- 
metry constraints (i) X is not a function of k y ; Y is not a 
function of k x (ii) X and Y are odd functions. Denoting 
lattice sites as j and j', 

H R = 

■^12 (j2 Yk v s[n ( k v(yj -vy)) s *$*>A«(' T *).s c f<>' 

- J2 Xk * ( kx ( X J ~ X 3')) S Vj >Vt> c )a {°v) „ a > c iV J • 

(110) 



Symmetry on a square lattice imply X u = Y u . The 
nearest-neighbor approximation involves keeping only 
terms with Xj — Xji = ±1 and yj — yy = ±1. In combina- 
tion, these are equivalent to setting X u = Y U = a R svo.u. 
Hence we arrive at Eq. (|26|) . 



B. Why minimize with respect to a and j3 

In Eq. (|47|) and (|48|) . we minimized G with respect to 
Ukri on the respective Fermi surfaces. This is equivalent 
to minimizing (H) — 7 (a 2 + 1 ) with respect to a and /3 
for the following reasons: 

(i) Keeping a 2 + /3 2 constant is equivalent to fixing 
particle number. This is strictly true in one dimension 
and approximately true in two dimensions for small to 
medium polarization. 

(ii) The Hohenberg-Kohn Theorem states that the 
ground state of a many-body Hamiltonian satisfies 

p{r) V(r) (111) 

where p is the density of electrons in real space and V 
is an external potential. Treating n^ri as a density in 
momentum space and Ek v (a, /3) as the external potential, 
a generalization of the Hohenberg-Kohn Theorem yields 

n kv <^ E kn <^=4> a,/3. (112) 



C. Critical J at small filling 

At J = J c , 

a 2 + P 2 = a 2 a (113) 

In the small filling approximation, the Fermi surfaces are 
circular, hence 

1 f^ F f 2n 
a Q ss— / dkk d8 cost kcosO) 
k 2 Jo Jo 

=- / dkkJ (k) 
Ti" Jo 
2 

— kp J\(Jvp^ 

TT 

=L y \ t) k 2m + 2 . (114) 

n /-^ 4 m (2m + 2 (m!) 2 F v ; 

m— 

J n (x) are Bessel functions of the first kind. 
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D. Connection with Reference |l4j 



E. Proof that a 2 + ■y 2 is a constant 



In Ref. [14[ the spin-split state was studied with 
a reduced Hamiltonian that limited interactions to 
only nearest- neighbor exchange (J) of electrons with 
antiparallel spin. We apply the formalism in this paper 
and derive an important result from Ref. fill ], namely 
the one-dimensional dispersion of the spin-split state. 
We derive this in the Rashba basis. 



We define operators h n = c^c,,, n = Y^n^v an< ^ 
A = Yliri r l? L ri- Following Ref. [T3 |. we only allow in- 
teractions in H { j =0) (Eq. (|19p ) between electrons with 
anti-parallel spin. We may deduce from Fig. (7) that 
interactions between electrons with antiparallel spin are 
equivalent to a sum of (i) interactions between electrons 
of the same Rashba polarization if both electrons have 
momenta with the same sign (ii) interactions between 
electrons of opposite Rashba polarizations if both elec- 
trons have momenta with opposite signs. We change Eq. 
(|19[) by making the substitution 

. . 1 + sgnO)sgn(fc') . „ „ . 

n k >n k — > [n k > + n k+ +n k >_n k -) 



1 — sgn(fc)sgn(fc') 



1 „ A sgn(fc)sgn(fc') 
---jn k >n k H 



(hk'+hk- + h k i_h k+ ) 

Afe'Afc. (115) 



After mean-field decoupling, the single-particle energy 
spectrum is 

(k) — — 2 — — sin k F cos e ^ cos k 

+ 2rj ( — sinfcj-sine jsin|fe|. (H6) 



By comparing Eq. (|36p with Eq. (Ill 6[) . we understand 
the spin-split Hamiltonian studied in Ref. 14| as the 
particular choices V = J and At = in the more general 
reduced Hamiltonian considered here. 



Consider an infinitesimal variation of a (defined in Eq. 
((BHI)) due to a relative displacement of the Fermi surfaces 



5a = — cos fc„ 5 



n ka 



(117) 



Since the shapes of the Fermi surfaces are unchanged for 
both spin up and down, this is mathematically equivalent 
to keeping the Fermi surfaces stationary and instead dis- 
placing the origin of the Brillouin zone in opposite direc- 
tions for spin up and down. We parametrize this displace- 
ment of the origin for spin {a = ±1} as 5k a — at (x+y). 
We may rewrite Eq. (|117[) as 

6a = — (cos {k v + ere) - cos k v J n ka . (118) 

To first order in e 

5a = —— ^2 sinfc^ a n ka = -e 7. (119) 



From Eq. (J67J), an analogous treatment of £7 reveals 

5j = ea. (120) 

Together, 



5a 7 
c>7 a 



(121) 
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